clear; clc; close all;
% 给出了 NMSE 与初始条件 z0 的关系

% 参数设置
mu = 3.5;
k = 0.16;
x0 = 0.1;
y0 = 0;

z0_range = linspace(-10, 8, 1000);

% 迭代步数设置
N_transient = 1000;   % 瞬态迭代步数
N_compute   = 1000;   % NMSE 计算步数

% 初始化存储 NMSE 数值
nmse_values = zeros(size(z0_range));

% 对不同的 z0 值调用 nmse.m 计算 NMSE
for i = 1:length(z0_range)
    z0 = z0_range(i);
    nmse_values(i) = nmse(x0, y0, z0, mu, k, N_transient, N_compute);
end

% 绘制 NMSE 曲线
figure;
plot(z0_range, nmse_values, 'LineWidth', 1.5);
xlabel('$z_0$', 'Interpreter', 'latex', 'FontSize', 14);
ylabel('NMSE', 'FontSize', 14);
title('NMSE 与 z_0 的关系');
grid on;
